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Abstract. Recent measurements on high-order harmonic generation (HHG) from N2 
aligned perpendicular to the driving laser polarization [B. K. McFarland el al, Science 
322, 1232 (2008)] have shown a maximum at the rotational half-revival. This has 
been interpreted as the signature of the HHG contribution from the molecular orbital 
just below the highest occupied molecular orbital (HOMO). By using the recently 
developed quantitative rescattering theory combined with accurate photoionization 
transition dipoles, we show that the maximum at the rotational half-revival is indeed 
associated with the HOMO-1 contribution. Our results also show that the HOMO-1 
contribution becomes increasingly more important near the HHG cutoff and therefore 
depends on the laser intensity. 
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High-order harmonic generation (HHG) has been extensively investigated both 
experimentally and theoretically over the last two decades [1]. Until very recently, 
the HHG has been understood as being due to tunneling ionization of an electron from 
the highest occupied molecular orbital (HOMO) and recombining back to the HOMO. 
The contribution from lower molecular orbitals are routinely neglected. That is not 
surprising since tunneling is a highly nonlinear process, and therefore highly selective to 
the HOMO due to energy considerations. In general, the neglect of the contribution from 
lower molecular orbitals are not justified for systems where the HOMO and HOMO-1 
are nearly degenerate, in other words, when the energy gap between the HOMO and 
lower molecular orbitals are much smaller than the ionization potential from the HOMO. 
Furthermore, for some molecular alignments, tunneling ionization from the HOMO is 
suppressed due to symmetry of the wavefunction [2] . Clearly, in that case the neglect of 
lower molecular orbitals is questionable. These two favorable conditions for observing a 
HOMO-1 contribution are present in N2, where the Itt^ HOMO-1 has a binding energy 
of 16.93 eV, quite close to the binding energy of the Sag HOMO (15.58 eV). 

Early theoretical calculations based on the strong-field approximation (SEA) model 
[3], H] have shown that harmonic yields from aligned N2 are maximum if the molecules are 
aligned along the laser polarization direction. These results are in good agreements with 
the pump-probe delay time experimental data [5l [6] as well as the recent more direct 
measurements [7j. On the other hand, the contribution from the HOMO-1 is expected 
to peak near 90°. Although the two molecular orbitals contribute to different alignment 
regions in the total harmonic yields, it is still a very challenging task to disentangle 
the HOMO-1 since its contribution is expected to be relatively weak. This is in strong 
contrast to the traditional single-photon photoionization where electrons are generally 
ionized from many MOs with comparable strengths (see, for example, [8]). 

In a recent experiment McEarland et al [9] reported that they have successfully 
observed the contribution from the HOMO-1 in aligned N2. That has been achieved 
within the pump-probe scheme with perpendicular pump-probe polarizations. Eor 
low harmonic orders below H23, the harmonic signals behave similarly to inverted 
< cos"^ 6 >, i.e., inverse of the degree of molecular alignment or 1— < cos^ 6* >. Eor 
higher harmonics, McEarland et al observed a maximum at the rotational half-revival, 
where inverted < cos^ 6* > is minimum. Eurthermore, the maximum at the half-revival 
is found to be quite pronounced in the HHG cutoff region, the location of which depends 
on the intensity of the driving laser. 

The goal of this paper is to show theoretically that the main features observed by 
McEarland et al [9] are indeed the signature of the HOMO-1 contribution. To support 
our claim, we have carried out calculations by using the recently developed quantitative 
rescattering theory (QRS) [10]. The photoionization transition dipole and its phase are 
obtained from state-of-the-art molecular photoionization calculations [H]. The QRS 
theory is based on the rescattering picture and it has been shown to give accurate 
results comparable with that from the time-dependent Schrodinger equation (TDSE) 
for rare-gas atoms [HI [13] and the molecular ion H^ [15]. The QRS has also been 
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Figure 1. Ionization rates from HOMO, HOMO-1, and HOMO-2 at laser intensities 
of 1 X 10" (a) and 2 x 10^^ W/cm^ (b). The data from HOMO-1 and HOMO-2 have 
been multipUed by a factor of 5 and 20, respectively. Calculations are carried out 
within the MO-ADK theory. 

shown to be able to reproduce most of the available experiments on aligned molecules 
CO2, O2 and N2 [inilin]- Other applications of the QRS include calculations of high- 
energy above-threshold ionization momentum and energy spectra (see, for example, [17] 
and references therein) and nonsequential double ionization [18] . Analytical derivations 
for the QRS have been reported quite recently in Refs. [191 1201 E]. In this paper we 
extend the QRS theory [10] to the multi-channel case, where the contribution from each 
channel to the total HHG induced dipole are added coherently. We assume that the ion 
cores are frozen during the time interval between ionization and photo-recombination. 

First we note that the HOMO (Scr^), HOMO-1 (l7r„) and HOMO-2 (2a„) have 
binding energies of 15.58, 16.93 and 18.73 eV, respectively. In order to have an idea 
about the magnitude of the relative contributions from the HOMO, HOMO-1 and 
HOMO-2 in the returning electron wave packet we compare the ionization rates from 
these molecular orbitals. These comparisons are shown in Fig. 1 for two different laser 
intensities of 1 x 10^^ and 2 x 10^"* W/cm^. The calculations were performed within 
the molecular tunnelling (MO-ADK) theory [2]. The MO-ADK C/ coefficients were 
calculated from asymptotic wavefunctions of each orbital, which in turn were obtained 
from the GAUSSIAN code [22]. Clearly the ionization rates depend strongly on the 
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Figure 2. Differential pliotoionization cross sections vs alignment angle from HOMO 
(a) and HOMO-1 (b) for some selected energies (expressed in units of harmonic orders 
for an 800-nm laser). 



alignment and the alignment dependent rates are different for different symmetries of 
the orbitals. This fact has been known before for HOMOs both theoretically [23l [2^ 
and experimentally [251126]. For the HOMO-1 (Ivr^j) the ionization rate peaks near 90°. 
Within this range of laser intensity, the rate from the HOMO-1 is a factor of 5 smaller 
than that of the HOMO, even at 90°. For the HOMO-2 (2cr„) the ionization rate peaks 
near 0°, similar to that of the HOMO. However its magnitude is significantly smaller 
than that of the HOMO (approximately by a factor of 20 for this range of laser intensity). 
We note that the relative contributions from the lower orbitals become increasingly more 
important as the laser intensity increases. We also carried out calculations using the 
SFA and found similar angular dependence for ionization rates from these three MOs. 
The relative strengths are also in good agreements with the MO-ADK, although the 
rate for the HOMO-2 from the SFA seems to be enhanced by a factor of 5, (i.e., still 
about a factor of 10 weaker than the HOMO). We note that our results do not agree 
with the recent results from the time-dependent density-functional theory calculations 
|27] . which show a peak in ionization rate from the HOMO-1 at an intermediate angle. 

The above analysis shows that the contribution to the returning wave packet from 
the HOMO-1 is about a factor of 5 smaller than that of the HOMO even near its 
peak at 90°. Can the HOMO-1 contribution to the harmonic generation be comparable 
with that of the HOMO? To answer this question we show in Fig. 2 a comparison 
of the differential photoionization cross sections from the HOMO and HOMO-1. (In 
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Figure 3. Contributions to HHG yields from HOMO and HOMO-1 for some selected 
harmonics with the laser intensity of 1.5 x 10^^ W/cm^ (left column) and 2 x 10^^ 
W/cm^ (right column). Molecular axis is assumed to be fixed (no averaging over 
alignment distribution was carried out). 



this paper we limit ourselves to the harmonics with polarization parallel to that of 
the driving laser. Therefore the relevant differential cross sections are for the electron 
emitted along the laser polarization direction [lO].) For convenience we express photon 
energy in units of the photon energy of the 800-nm laser (1.55 eV). First we note that the 
cross sections from the HOMO vary strongly from one energy to the next. Nevertheless, 
one general feature can be seen is that the cross sections at large angles near 90° are 
quite small, say, about 0.25 Mb for H29, and decreases quickly with energy. On the 
other hand the photoionization cross sections from the HOMO-1 all have a dominant 
peak at 90°, which reaches 1.25 Mb for H29. Recall that according to the QRS theory, 
HHG yield is proportional to a product of the returning wave packet and the differential 
photoionization cross section. The results shown in Figs. 1 and 2 therefore indicate that 
the HOMO and HOMO-1 contributions could be comparable near 90°. 

Having established qualitatively that the contribution from the HOMO-1 cannot 
be neglected for large angles, we now analyze the HHG yields from the actual QRS 
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Figure 4. Harmonic signals as functions of delay time near half-revival. The laser 
intensity is 2 x lO^** W/cm^ (a) and 1.5 x 10^'' W/cm^ (b). A (pump) laser pulse of 
90 fs duration (FWHM), intensity of 3 x 10^'^ W/cm^ is used to align the molecules. 

calculations. In Fig. 3(a-d) (left column) we show the HHG yields from the HOMO 
and HOMO-1 for H19, H29, H33 and H35. The 800-nm laser pulse is of 30fs duration 
(FWHM) and intensity of 1.5 x 10^^ W/cm^. Clearly, the HOMO contribution dominates 
for alignment angles smaller than 45° for all the harmonics. For large angles the HOMO- 

1 becomes comparable with the HOMO already near H25 and dominates for the higher 
harmonics, especially beyond the cutoff at H29. Similar pattern repeats at a higher 
laser intensity of 2 x 10^^ W/cm^, shown in Fig. 3(e-h) (right column). However the 
HOMO-1 is only comparable with the HOMO at large angles for H29 and it starts to 
dominate only at higher harmonics. The enhanced contribution from the HOMO-1 near 
the cutoff can be understood as the consequence of the delay in the harmonic cutoff for 
the HOMO-1, since the ionization potential from the HOMO-1 is about 1.5 eV greater 
than that from the HOMO. As the cutoff moves to near H35 for the laser intensity of 

2 X 10^^ W/cm^, the HOMO-1 contribution dominates at much higher harmonic orders 
as compared to the lower intensity case shown in the left column. This fact has been 
noticed earlier by McFarland ei a/ [9]. 

The main results of this paper are presented in Fig. 4, where we show the harmonic 
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signals as functions of delay time between the pump and probe laser pulses at two 
intensities of 1.5 x 10^^ and 2 x 10^^ W/cm^. The signals have been normalized to 
the signals from the isotropic distribution. These results can be compared directly 
with the measured HHG signals by McFarland et al Similarly to the experiments 
by McFarland et al j9], we take the pump laser polarization to be perpendicular to 
the probe laser polarization. This is used in order to facilitate the observation of the 
harmonic signals at the large angles where the HOMO-1 is dominant for high harmonic 
orders. Theoretically the harmonic signals are obtained by a coherent convolution of 
the HHG yields from the QRS calculations with the partial alignment distribution. The 
time-dependent molecular alignment distribution is calculated by solving the TDSE for 
N2 molecules in the pump (alignment) laser field within the rotor model |28]. The pump 
(alignment) laser has a pulse length of 90 fs (FWHM), intensity of 3 x 10^^ W/cm^ and 
800-nm wavelength. We assume the Boltzmann distribution for the rotational levels 
at the initial time and the rotational temperature is taken to be 40 K. The above 
parameters are chosen to closely match the experimental conditions of McFarland et al 

At lower laser intensity of 1.5 x 10^^ W/cm^, as can be seen from Fig. 4(b), the lower 
harmonic (as represented by H17 and H19) has a minimum at the half-revival near 4.1 
ps. In other words, the HHG signal behaves as inverted < cos^ 6* > (not shown), which 
measures the degree of molecular alignment. For H29 and higher harmonics a peak 
superimposed on the minimum can be seen. This is clear evidence for the increasing 
importance of the contribution from the HOMO-1 at large angles, shown in Fig. 3(a). 
Recall that at the half-revival near 4.1 ps, the molecules are maximally aligned along 
the pump polarization direction, which is perpendicular to probe polarization. If we 
artificially remove the HOMO-1 contribution, all the harmonics behave similarly to HI 7, 
i.e., as inverted < cos^ 9 >. On the other hand, the contribution from the HOMO-2 are 
found to be negligibly small and its inclusion does not affect our results. Similar behavior 
is seen at higher laser intensity of 2 x 10^^ W/cm^ shown in Fig. 4(a). However the peak 
at 4.1 ps starts to show up only at H31 and more systematically after H35. This is due to 
the fact that the harmonic cutoff is shifted to H35 at this high intensity. It is clear that 
our QRS theoretical results for a single-molecule response already reproduce quite well 
the general behavior reported by McFarland et a/ [9]. We also performed calculations 
using the standard SFA, with and without including the multiple orbitals, and found 
that the SFA does not reproduce experimental results. For a more complete theory, 
one certainly needs to include the macroscopic propagation. As has been reported quite 
recently by Sickmiller and Jones [29], the phase-matching condition can have a dramatic 
effect on the HHG signal in the waveguide. The effect of propagation is expected to be 
much less significant in the gas jet experiment by McFarland et al [9] than in Ref. |29j . 
At present the macroscopic propagation of the HHG in molecular gas media is still 
largely unexplored with only one attempt for H.2 [30] having been reported so far. 

In conclusion, by using the recently developed quantitative rescattering theory 
(QRS) combined with accurate photoionization transition dipoles we have confirmed 
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theoretically that the peaks superimposed on the minimum near rotational half- 
revival observed experimentally in N2 by McFarland et al [9] can be attributed to 
the contribution from the HOMO-1. Our results also show that the contribution 
from the HOMO-1 becomes more important in the HHG cutoff region and therefore 
depends on the laser intensity. In this paper we have limited ourselves to the parallel 
polarization component for the emitted harmonics, as only the case of orthogonally 
polarized pump and probe pulses was considered. A recent experiment by Zhou et al 
[3T] and an earlier experiment by Levesque et al [32], in which the angle between pump 
and probe polarizations were varied, showed that the ratio of the perpendicular to 
parallel component intensities could be as high as ~ 0.2. However this ratio decreases 
significantly for large angles between the pump and probe polarizations and goes to 
zero at 90° due to the symmetry. The importance of this effect is currently under 
investigation. Finally we mention the recent experimental and theoretical results by 
Smirnova et al [33], which presented evidence for the multiple orbitals effect in CO2. 

We thank M. Giihr, B. McFarland, J. Farrell, and P. Bucksbaum for communicating 
their experimental data to us and for valuable discussions. This work was supported 
in part by the Chemical Sciences, Geosciences and Biosciences Division, Office of Basic 
Energy Sciences, Office of Science, U. S. Department of Energy. 
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